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for  a specific  design  example. 
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NOMENCLATURE 


= propeller  thrust  coefficient 


,SPVa2ttR2 


fo/c 


propeller  diameter 

blade  section  camber/chord  ratio 


advance  coefficient  _A 

nD 


= number  of  propeller  blades 


torgue  coefficient 
thrust  coefficient 


pn2D5 


pn2D4 

propeller  rotational  speed  in  revolutions/second 

blade  section  pitch 

pitch  of  trailing  vortex  sheet 


pitch  of  trailing  vortex  sheet  according  to  lifting 
line  theory 
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U„rU4 


ua'ut 


V. 


pitch  of  trailing  vortex  in  ultimate  wake 

propeller  torque 
radial  coordinate 
tip  vortex  core  radius 

radium  of  curvature  of  helical  vortex 

radius  of  ultimate  wake  helix 

propeller  thrust 

axial,  tangential  induced  velocities 

normalized  axial  and  tangential  induced  velocities 

propeller  speed  of  advance 


geometric  advance  angle  tan 


-1 


2wnr 


advance  angle  including  induced  velocities  from 
lifting-linfe  theory 

advance  angle  of  vortex  ultimate  wake 
circulation 

advance  coefficient  g-tanB. 

K 3- 

K J 

propeller  efficiency  T 

2„Kq 

fluid  mass  density 

propeller  rotational  speed  in  radians/second 
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I. 


INTRODUCTION 


A major  obstacle  in  obtaining  an  "exact"  numerical  solution 
of  the  flow  around  a three-dimensional  lifting  surface  is  the 
difficulty  in  determining  the  position  of  the  trailing  vortex 
wake.  This  problem  has  the  usual  difficulties  associated 
with  boundary  value  problems  with  unknown  boundaries,  with  the 
additional  complication  of  the  essential  role  of  viscosity 
in  the  tip  vortex  core  region. 

It  is  fortunate  that  for  many  applications  a precise 
knowledge  of  the  geometry  of  the  trailing  vortex  wake  is  not 
essential.  The  pressure  distribution  on  a three-dimensional 
planar  hydrofoil  can  be  accurately  computed,  in  most  cases, 
with  complete  neglect  of  the  deformation  of  the  wake.  Exceptions 
to  this  statement  include  the  determination  of  the  pressure 
distribution  in  the  immediate  vicinity  of  the  tip,  or  the 
analysis  of  low  aspect  ratio  hydrofoils  operating  with  vortex 
sheet  separation  from  the  leading  edge. 

The  problem  is  more  complicated  in  the  case  of  a propeller. 
Here  the  trailing  vortices  follow  a generally  helical  path 
downstream,  and  the  velocity  which  they  induce  on  the  blade 
is  significantly  altered  by  a deformation  in  their  trajectory. 

For  example,  the  velocity  induced  on  one  blade  by  the  tip 
vortex  shed  from  the  next  blade  ahead  is  roughly  proportional 
to  the  inverse  of  the  minimum  distance  between  the  blade  and 
the  vortex.  This  distance  is  significantly  altered  by  the 
combination  of  axial  and  radial  deformation  of  the  vortex. 

This  problem  had  been  recognized  in  the  earliest  analytical 
treatments  of  propeller  theory,  and  the  well-known  represen- 
tation of  the  "moderately  loaded  propeller"  can  be  thought  of 
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as  a deformed  wake  model  [1],  12].  In  this  case,  a very 
simple  deformation  is  introduced  by  modifying  the  pitch  of 
the  trailing  vortices  to  match  the  flow  angle  computed  at  a 
lifting-line  with  the  same  radial  circulation  distribution 
as  the  actual  propeller.  Since  the  trailing  vortices  shed 
at  different  radii  will  not  necessarily  have  the  same  pitch 
the  resulting  vortex  sheet  for  downstream  will  be  of  con- 
tinuously changing  form. 

This  model  has  had  the  enviable  characteristic  of  getting 
better  with  age.  When  first  introduced,  computer-aided  numerical 
lifting-surface  theory  was,  of  course,  not  available  to  determine 
detailed  propeller  blade  section  characteristics  for  a prescribed 
loading.  Semi-empirical  camber  and  pitch  correction  factors 
were,  of  necessity,  employed  in  determining  blade  section 
shape  [3],  so  that  errors  in  the  final  result  could  not  necessarily 
be  blamed  on  the  wake  model.  As  numerical  lifting  surface 
theory  progressed  in  accuracy  [4] , [5] , [6] , [7] , correlation 
between  theory  and  experiment  became  excellent,  even  though 
the  classic  "moderately-loaded"  wake  model  was  retained. 

It  is  important  to  recognize,  however,  that  this  model 
cannot  properly  represent  flow  conditions  downstream  of  the 
propeller.  The  author  must  admit  falling  into  this  trap  in 
developing  a numerical  procedure  for  computing  propeller  field 
point  velocities  [8],  [9],  Retention  of  the  "moderately- 
loaded"  wake  model,  in  this  case,  results  in  an  oscillation  in 
the  amplitude  of  the  blade  freqeuncy  pressure  with  downstream 
positions,  while  both  experiment  and  pure  linear  theory  yield 
monotonically  decreasing  amplitudes.  On  the  other  hand, 
correlation  between  theory  and  experiment  near  the  propeller, 
which  is  generally  of  primary  importance,  is  extremely  satisfactory. 

One  possible  variation  on  the  classical  "moderately-loaded" 
wake  model  was  introduced  by  Burrill[10].  Recognizing  that 
the  velocity  induced  by  a helical  vortex  is  doubled  downstream. 
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the  pitch  of  the  helical  vortices  shed  from  the  blades  was 
appropriately  increased.  Since  the  Bur  rill  procedure , which 
again  pre-dated  the  introduction  of  digital  computers,  was 
handicapped  by  possible  inaccuracies  in  estimating  lifting- 
surface  effects,  the  final  result,  as  with  [3],  could  not  be 
used  to  judge  the  effectiveness  of  the  wake  model. 

One  can  conclude,  therefore,  that  the  "moderately-loaded" 
wake  model  is  evidently  an  effective  "substitute  vortex  system" 
for  the  purpose  of  determining  the  flow  on  or  near  the  blade. 
Viewed  in  this  way,  it  does  not  matter  whether  or  not  the 
substitute  vortex  system  deviates  significantly  from  the  true 
vortex  wake  downstream.  Conversely,  since  it  is  a substitute, 
it  may  lead  to  wrong  answers  if  used  for  any  other  purpose. 

Considerable  effort  has  been  directed  in  recent  years 
towards  the  determination  of  the  true  propeller  wake  geometry. 
Since  a concurrent  survey  paper  by  Cummings  [11]  reviews  this 
work,  no  attempt  will  be  made  here  to  duplicate  this  material, 
except  to  state  the  general  conclusion.  Detailed  numerical 
computations  of  propeller  vortex  sheet  deformations,  while 
within  the  realm  of  possibility  for  current  computers,  are 
probably  not  worth  the  effort.  At  the  present  time,  we  are 
unaware  of  the  existence  of  any  practical  working  programs  for 
this  purpose. 

We  have  therefore  pursued  a simpler  approach  which  is 
intermediate  between  the  theory  of  moderately-loaded  propellers 
and  a complete  numerical  determination  of  the  wake  geometry. 
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II. 


THE  ULTIMATE  WAKE  MODEL 


Observations  of  propeller  tip  vortex  trajectories  in  the 
MIT  Variable  Pressure  Water  Tunnel  tend  to  indicate  that  the 
roll-up  process  must  be  extremely  rapid.  The  slipstream 
contracts  behind  the  propeller  to  a constant  value  which  is 
maintained  for  the  entire  length  of  the  test  section  downstream 
of  the  propeller.  If  any  significant  amount  of  vortex  roll-up 
were  taking  place  further  downstream,  this  would  be  evidenced 
by  a continual  change  in  radius  and  pitch.  It  seemed,  therefore, 
that  rather  than  chasing  vortices,  it  would  be  more  efficient  to 
focus  attention  on  the  ultimate  wake.  The  ultimate  wake  would  consist 
of  a concentrated  helical  tip  vortex  shed  from  each  blade,  and 
a hub  vortex  coincident  with  the  shaft  axis.  Since  we  are 
looking  far  downstream,  these  vortices  can  be  regarded  as  infinite 
in  extent,  so  that  the  problem  becomes  effectively  a two- 
dimensional  one. 

The  ultimate  wake  must  be  joined  to  the  propeller  by  a 
"transition  wake"  region,  which  is  where  the  complex  rolling-up 
and  contracting  phenomenon  takes  pla^e.  However,  if  this 
region  is  limited  in  extent,  it  can  be  modelled  with  sufficient 
accuracy  by  a smooth  surface  consisting  of  vortex  lines  shed 
from  the  trailing  edge  and  converging  at  the  "roll-up  points". 

This  concept  is  illustrated  in  Figure  1. 

The  development  of  this  ultimate  wake  model  was  undertaken 
by  Loukakis  [12] , who  developed  analytical  methods  to  determine 
the  ultimate  wake  radius  and  pitch,  and  ran  a pilot  experiment 
in  the  water  tunnel  to  measure  these  quantities.  Good  agreement 
was  found  between  theory  and  experiment.  His  results  showed, 
in  particular,  that  the  ultimate  wake  radius  was  practically 
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Schematic  of  simplified  deformed  wake' model  consisting  of  a transition  and 
an  ultimate  wake. 


pitch  was  less  than  that  predicted  by  the  theory  of  moderately- 
loaded  propellers.  The  latter  tended  to  confirm,  at  least 
qualitatively,  the  similar  conclusion  reached  by  Cummings  [13]. 

The  present  research  is  directed  towards  the  application 
of  this  deformed  wake  model  to  design,  and  consists  of  three 
parts: 


1. 


A systematic  series  of  measurements  of  ultimate 
wake  pitch  and  radius. 


2. 

3. 


Correlation  of  these  measurements  with  theory. 
Modification  of  the  existing  propeller  lifting- 


surface  program  to  include  the  computation  of 
induced  velocities  on  the  blade  in  accordance 
with  the  present  deformed  wake  model. 


III.  EXPERIMENTAL  PROGRAM 


I: 


i 
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A.  Performance  Tests 

A series  of  four  three-bladed  propellers  were  supplied 
by  DTNSRDC  for  the  purpose  of  this  study.  The  characteristics 
of  the  propellers,  which  all  have  a pitch/diameter  ratio  of 
approximately  1.07,  are  given  in  [14],  Figure  2,  reproduced 
from  [14] , shows  the  variations  in  blade  outline  covered  by 
the  series. 

In  addition,  a stock  three-bladed  motorboat  propeller  with 
a pitch/diameter  ratio  of  0.67  was  tested  in  order  to  obtain 
data  applicable  to  propellers  with  a lower  design  advance 
coefficient. 

Measurements  of  thrust  coefficient,  Rj.  , torque  coefficient, 
Kq  , and  efficiency  versus  advance  coefficient  J , were  made 
in  the  MIT  Variable  Pressure  Water  Tunnel  under  non-cavitating 
conditions.  Following  standard  procedures,  the  test  RPM  was 
adjusted  to  provide  the  highest  Reynolds  Number  possible  without 
exceeding  the  capacity  of  the  dynamometer  or  the  maximum  tunnel 
flow  speed.  For  these  propellers,  the  range  of  RPM  was  from 
1000  to  1800  , yielding  Reynolds  Numbers  ranging  from  2 * 106 
to  3.5  x 10G. 

Faired  values  of  Kip  , Kq  , and  efficiency  corrected  for 
tunnel  wall  effect  at  intervals  of  0.1  in  advance  coefficient 
are  given  in  Tables  1-5.  The  results  for  the  four  series 
propellers  are  generally  in  good  agreement  with  the  results  of 
the  DTNSRDC  open  water  tests  [14]. 
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TABIC  1.  FAIRED  EXPERIMENTAL  VALUES  OF  PROPELLER  PERFORMANCE 
CORRECTED  FOR  TUNNEL  WALL  EFFECT  FOR  DTNSRDC  PROPELLER  4118. 


3 

•4 

n 

.420 

.0675 

.198 

.3 

.381 

.0618 

.295 

.4 

.342 

.0360 

.589 

.5 

.301 

.0301 

.478 

.6 

.237 

•0460 

.559 

• J 

.213 

.0377 

.828 

.8 

.170 

.0318 

.682 

.9 

.129 

.0260 

.714 

1.0 

.086 

.0193 

.700 

1.1 

.040 

.0126 

.558 

1.2 

-.005 

.0055 

-.224 

TABU  2.  FAIRED  EXPERIMENTAL  VALUES  OF  PROPELLER  PERFORMANCE 
CORRECTED  FOR  TUNNEL  WALL  EFFECT  FOR  DTNSRDC  PROPELLER  4112. 


3 

s 

*Q 

0 

.1 

.332 

•0472 

.227 

.3 

.310 

.0444 

.333 

.4 

.293 

.0416 

.434 

.5 

.255 

.0386 

.326 

.4 

.226 

.0333 

.608 

.7 

.195 

.0321 

.676 

.8 

.163 

.0283 

.729 

.9 

.133 

.0246 

.773 

1.0 

.100 

.0199 

.804 

1.1 

.064 

.0142 

.793 

1.2 

.024 

.0080 

.628 

1.3 

-.013 

.0014 

-1.904 

TABU  J.  FAIRED  EXPERIMENTAL  VALUES  OF  PROPELLER  PERFORMANCE 
CORRECTED  FOR  TUNNEL  WALL  EFFECT  FOR  DTNSRDC  PROPELLER  4133 


j 

*T 

n 

.2 

.523 

.0894 

.188 

.) 

.467 

.0803 

.277 

.4 

.411 

.0717 

.365 

• 3 

.334 

.0628 

.649 

.6 

.296 

.0336 

.324 

.2 

.237 

.0445 

.593 

.8 

.182 

.0341 

.645 

.9 

.130 

.0282 

.440 

1.0 

.079 

.0208 

.405 

1.1 

.028 

.0135 

.357 

1.2 

-.024 

.0042 

•.730 
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TABLE  4.  FAIRED  EXPERIMENTAL  VALUES  OF  PROPELLER  PERFORMANCE 
CORRECTED  FOR  TUNNEL  WALL  EFFECT  FOR  DTNSRDC  PROPELLER  4143. 


J 

*T 

n 

.419 

.0629 

.212 

.) 

.171 

.0566 

.314 

.4 

. 328 

.0506 

.412 

.292 

.0445 

.303 

.234 

.0383 

.sat 

.7 

.190 

.0323 

.434 

.a 

.148 

.0267 

.704 

.9 

.110 

.0217 

.724 

1.0 

.071 

.0167 

.482 

1.1 

.012 

.0116 

.487 

1.2 

•.007 

.0065 

.214 

TABLE  5.  FAIRED  EXPERIMENTAL  VALUES  OF  PROPELLER  PERFORMANCE 
CORRECTED  FOR  TUNNEL  WALL  EFFECT  FOR  12  x 8 STOCK  PROPELLER. 


J 


s 

*0 

n 

,244 

.0294 

.144 

,229 

.0244 

.261 

,192 

.0224 

.406 

,132 

.0190 

.510 

.113 

.0134 

.582 

,074 

.0114 

•605 

.030 

.0072 

.462 

.014 

.0023 

-.839 
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well  suited  for  the  measurement  of  propeller  tip  vortices, 
due  to  its  square  test  section,  with  one  long  unobstructed 
plexiglass  window  on  each  side.  The  propeller  can  be  located 
in  the  forward  portion  of  the  test  section,  thus  leaving 

sufficient  distance  downstream  for  observation  of  the  tip  j 

vortices.  The  tunnel  pressure  can  be  adjusted  to  provide 

the  minimum  amount  of  tip  vortex  cavitation  to  permit  observation. 

The  initial  experiment  conducted  by  Loukakis  [12]  was  done,  j 

of  necessity,  in  the  simplest  possible  way.  A sighting  device 
was  made  by  gluing  together  three  plastic  drafting  triangles. 

This  device  was  held  firmly  against  the  test  section  window, 
and  the  tip  vortices  were  traced  by  sighting  along  the  edge  of 

the  device  normal  to  the  window.  ’ 

j 

For  the  present  study,  an  improved  sighting  instrument 
was  devised.  This  consisted  of  a ten-power  telescopic  sight  j 

attached  to  a drafting  machine  which,  in  turn,  was  securely  \ 

attached  to  the  tunnel  test  section.  The  telescopic  sight 
could  be  moved  along  a tip  vortex  with  great  precision,  and 
the  horizontal  and  vertical  coordinates  of  a set  of  points 
could  be  read  directly  on  the  drafting  machine  scales.  A 
photograph  of  this  equipment  in  operation  appears  in  Figure  2. 

With  the  telescope  centered  on  a vortex,  it  was  immediately 

apparent  that  tunnel  pressure  could  be  changed  considerably  ! 

before  any  alteration  in  the  position  of  the  vortex  occurred.  j 

The  data  was  analysed  by  fitting  the  best  least-squares  ! 

helix  through  the  measured  points.  This  is  equivalent  to  the  ] 

problem  of  obtaining  the  amplitude  and  phase  of  a harmonic 
function  of  unknown  wave  length.  The  procedure  is  to  assume 
the  wavelength,  i.e.  the  pitch  of  the  helix,  solve  for  the 
amplitude  and  phase  of  the  fundamental  by  the  usual  technique 
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Figure  2.  - Photograph  of  telescopic  sighting  device  for  the 
measurement  of  tip  vortex  trajectories. 
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of  discrete  harmonic  analysis,  and  compute  the  mean-squared 
error.  The  wave-length  is  then  altered  slightly,  and  the 
harmonic  analysis  and  error  calculation  is  repeated.  An 
iterative  procedure  is  then  followed  to  determine  the  wavelength 
for  minimum  error. 

A typical  computer  output  is  shown  in  Figure  3.  This 
consists^of  a tabulation  of  input  coordinates,  a plot  of  the 
original  data  together  with  the  least-squares  helix,  and  a 
tabulation  of  the  final  pitch,  radius  and  phase  of  the  helix. 

In  almost  all  cases,  the  helix  provided  an  extremely  close 
approximation  to  the  measured  data. 

It  was  found  that  the  accuracy  of  these  measurements  was 
very  dependent  upon  the  quality  of  the  propeller  models.  If 
the  blades  were  not  all  identical,  the  wake  would  become  progres- 
sively more  distorted  downstream.  The  hub  vortex,  in  this  case, 
would  not  remain  on  the  axis,  but  would  be  attracted  towards 
one  tip  vortex,  and  in  the  process  would  develop  a helical  shape 
with  increasing  radius.  Eventually,  this  would  cause  complete 
breakdown  of  the  original  tip  vortices. 

While  some  laboratory  model  propellers  tested  were  far 
more  symmetrical,  the  differences  in  pitch  of  the  tip  vortices 
from  each  blade  could  readily  be  measured.  Some  models  were 
clearly  superior  to  others  in  this  respect.  This  effect  was 
particularly  noticeable  for  the  stock  motorboat  propeller,  where 
one  would  not  expect  laboratory  standard  tolerances. 

To  obtain  the  best  average  value  of  the  ultimate  wake  pitch 
and  radius,  measurements  were  made  of  the  tip  vortex  from  each 
blade  set  at  both  zero  and  ninety  degrees  by  means  of  the 
phase  of  the  strobe-light. 

The  results  are  given  in  Figures  4-7  for  a series  of  four 
three-bladed  propellers  designed  by  DTNSRDC  [14],  The  curve 
labeled  "initial  pitch"  is  the  axial  distance  from  the  tip  of 
the  blade  to  the  first  revolution  of  the  tip  vortex.  The 
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Figure  3.  - Sample  computer  output  of  vortex  curve-fitting  procedure. 


ultimate  pitch  Pw  , is  the  pitch  of  the  tip  vortex  measured 


as  far  downstream  in  the  test  section  as  possible.  At  the 
design  advance  coefficient,  J=  0.833  , the  average  initial 
pitch  was  three  percent  less  than  the  ultimate  pitch.  This 
extremely  small  difference  is  one  indication  that  vortex  sheet 
roll-up  is  extremely  rapid,  and  that  one  does  not  have  to  proceed 
very  far  downstream  to  reach  the  ultimate  wake. 

The  two  dotted  lines  plotted  in  Figures  4-7  are,  respectively, 
the  pitch  of  the  undisturbed  flow 

(1) 


PQ  = 2irr  tan  B , 

and  the  pitch  of  the  tip  vortex  in  accordance  with  the  theory 
of  moderately  loaded  propellers 

P^  = 2irr  tan  B^  . (2) 


The  latter  was  obtained  from  a lifting-line  theory  calculation 
with  the  measured  thrust  as  input,  assuming  that  the  radial 
distribution  of  circulation  was  optimum. 

In  all  cases,  the  ultimate  wake  pitch  falls  almost  exactly 
midway  between  these  two  limits. 

The  ultimate  wake  radius  is  plotted  at  the  bottom  of 
Figures  4-7  and  can  be  seen  to  be  essentially  independent  of 
propeller  loading,  with  a mean  value  of  83%  of  the  propeller 
radius.  Subsequent  observations  of  a number  of  propellers  with 
varying  geometry  and  number  of  blades  have  been  made,  and  the 
ultimate  wake  radius  is  generally  very  close  to  83%  . However, 
at  extremely  high  trust  coefficients,  one  can  observe  a dependence 
of  ultimate  wake  radius  on  loading,  with  values  falling  below 
83%  . Some  indication  of  this  trend  is  apparent  from  Figures 
4-7. 
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IV.  COMPUTATION  OF  THE  ULTIMATE  WAKE  PITCH 

A.  Discrete  Helical  Vortices 

It  is  assumed  that  all  vorticity  rolls  up  into  discrete 
tip  vortices  and  a hub  vortex,  in  a relatively  short  distance 
behind  the  propeller.  In  the  ultimate  wake,  we  then  have  a 
flow  consisting  of  set  of  infinite  helical  vortex  cores,  and 
an  infinite  straight  hub  vortex  core.  The  core  radius,  rc  , 
is  assumed  to  be  known. 

The  self-induced  velocity  of  a set  of  K infinite  helical 
vortex  can  be  expressed  as  follows, 


u 

a 


Ua(ew'rc'K) 


ut  = ut(Bw'  rc'  K) 


(3) 


where  is  the  radius  of  the  helix,  with  pitch 

Pw  = 2 7i tan  $w  . (4) 

r is  the  circulation,  and  rc  is  the  core  radius  non-dimension- 

alized  on  the  radius  of  the  helix.  The  functions  u and  u. 

a t 

contain  a logarithmic  singularity  as  rc  -*■  0 , so  that  the 
assumption  of  a finite  core  size  is  essential. 

Numerical  values  for  and  ufc  may  be  computed  in 

two  parts.  The  part  of  the  helix  close  to  the  point  where  we 
are  computing  the  velocity  is  approximated  by  a segment  of  a 
circular  vortex  with  radius  equal  to  the  radius  of  curvature 
of  the  helix.  As  shown,  for  example,  by  Lamb  [15] , this  yields 
a velocity  normal  to  the  helix  of  magnitude 
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6u  = 
r 


where  2<J>  is  the  included  angle  of  the  segment,  and  RR  is 
the  radius  of  the  tangent  ring, 

RR  = Vcos2f3w  * (6) 

The  velocity  induced  by  the  rest  of  the  first  helix,  and  from 
the  entire  remaining  K-l  helices,  may  be  obtained  by  numerical 
integration  of  the  usual  Biot-Savart  Law  integral  for  zero  core 
size.  Numerical  experimentation  indicates  that  the  answer  is 
insensitive  to  the  choice  of  2<j>  when  the  length  of  the  ring 
segment  is  varied  from  4-20  times  the  core  radius. 


Calculations  of  u_  and  u. 


for  various  values  of  Bw  , 


r and  K are  given  by  Loukakis  [12]  in  Tables  6-7  and  Figures 

C 

8-9,  which  are  reproduced  here. 

The  hub  vortex  induces  a tangential  velocity,  independent 
of  core  size 

u = ~K^  (7) 

t 2ttRw  * Ui 


The  requirement  that  the  tip  vortex  helix  be  a streamline 
leads  to  the  following  equation  for  8W  : 


tan  8W  = 


VA  + 5a(BW'  rc'  K>  • r/*w 

-KF  — 

“Rw-  ns;  + VSw'rc'K>  • 


The  ultimate  wake  pitch  may  therefore  be  computed  by  an 
iterative  procedure  from  (8)  once  the  wake  radius  and  core 
size  is  given.  Loukakis  [12]  has  shown  that  if  the  propeller 
thrust  is  known,  a momentum  balance  can  be  used  to  obtain  one 
more  relationship  between  pitch  and  radius,  thus  permitting  a 
solution  for  both  quantities.  However,  we  can  avoid  this 
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Table  6.  - Self-induced  velocities  as  a function  of  number  of  blades.  Ratio  of  core 
radius  to  vortex  radius  is  1/50.  [From  Loukakis  (12)] 
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Table  7.  - Self-induced  velocities  as  a function  of  core  radius.  [From  Loukakis  (12)] 


RATIO  = 


complication  if  we  are  willing  to  assume  a value  for  ultimate 
wake  radius  at  the  outset.  On  the  basis  of  our  experiments, 
it  would  seem  that  this  could  be  assumed  to  be  83%  of  the 
propeller  radius  for  all  applications  to  light  or  moderately 
loaded  propellers. 

In  either  case,  we  must  still  assume  a core  radius.  The 
sensitivity  of  this  variable  can  be  judged  from  Loukakis' 
calculations  plotted  in  Figures  8-9  for  assumed  values  of  core 
radius  ranging  from  one  to  ten  percent  of  the  vortex  radius. 

Since  the  contribution  of  only  one  blade  is  affected,  precise 
knowledge  of  the  actual  core  radius  is  not  critical,  particularly 
for  propellers  with  large  numbers  of  blades. 

To  compare  results  obtained  in  this  way  with  the  experiments, 
it  is  necessary  to  make  allowance  for  the  influence  of  the  tunnel 
walls  on  the  pitch  and  radius  of  the  ultimate  wake. 

An  estimate  of  the  change  in  pitch  can  be  made  by  considering 
the  backflow  velocity  induced  by  the  tunnel  walls.  This  backflow 
evaluated  in  the  plane  of  the  propeller  is  the  basis  for  cor- 
recting the  advance  coefficient  in  the  tunnel  to  match  that  of 
an  open-water  test.  This  correction  was  derived  for  an  actuator 
disc  by  Glauert  [16].  B.  D.  Cox  [17]  computed  the  backflow 
for  a lifting-surface  representation  of  the  propeller  using  the 
MIT  field  point  program,  and  found  that  for  the  large  clearances 
present  in  a usual  tunnel  test  that  the  results  were  essentially 
the  same.  The  classical  Glauert  correction  has  therefore  been 
built-in  to  the  propeller  data  reduction  program  used  at  MIT. 

The  corrected  tunnel  advance  coefficient  equalizes  the 
inflow  in  the  plane  of  the  propeller  in  the  tunnel  and  in 
open  water.  However,  since  the  backflow  velocity  doubles 
downstream,  the  pitch  of  the  ultimate  wake  will  be  reduced  by 
the  ratio  of  this  change  in  backflow  to  the  velocity  in  the 
plane  of  the  propeller.  This  correction,  which  is  a function 
of  propeller-loading,  may  be  obtained  from  the  processed 
performance  data. 
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An  estimate  of  the  influence  of  the  tunnel  walls  on  the 
radius  of  the  ultimate  wake  can  also  be  made  by  comparing  the 
results  of  an  actuator  disc  in  an  unlimited  fluid  with  an 
actuator  disc  in  an  infinite  cylinder.  The  following  results 
are  obtained  for  a 12-inch  diamter  propeller  in  a 20x20 
inch  test  section  with  a thrust  coefficient  CT  = 0.565  . 

Actuator  Disc  in  Unlimited  Fluid  Ry^/R  = 0.95 

" " in  a Tunnel  " =0.96 

Observed  in  Tunnel  " =0.83 

While  the  tunnel  walls  inhibit  slipstream  contraction  slightly, 
this  influence  is  negligible  in  comparison  with  the  influence 
of  vortex  roll-up,  and  is  therefore  not  considered  further. 


§- 


B.  Ultimate  Wake  Pitch  Determined  from  Actuator  Disc  Theory 

While  the  actuator  disc  does  not  provide  a suitable 
estimate  of  ultimate  wake  radius  due  to  vortex  roll-up,  it  can 
provide  an  extremely  simple  estimate  of  the  pitch.  We  can 
consider  the  actuator  disc  to  be  represented  by  a tube  of  ring 
vortices  proceeding  downstream  from  the  plane  of  the  disc. 

It  can  be  shown  from  energy  considerations  that  this  vorticity 
must  be  convected  at  the  average  of  the  slipstream  velocity 
and  the  free-stream  velocity. 

If  we  assume  that  this  convection  velocity  is  the  same 
for  a propeller  with  a finite  number  of  blades  operating  at 
an  advance  coefficient,  J , we  obtain  the  simple  expression 
for  the  pitch  of  the  ultimate  wake 
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(9) 
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V. 


RESULTS 


A summary  of  measured  and  calculated  results  for  the 
DTNSRCD  three-bladed  propeller  series  is  given  in  Table  8. 

Two  advance  coefficients  were  considered;  the  design  value  of 
J = 0.833,  and  a lower  value  of  0.6.  The  latter  results  in  a 
thrust  coefficient,  CT  , of  about  three  times  the  design  value. 

At  the  low  advance  coefficient,  the  simple  actuator  disc 
formula  (10)  predicts  values  of  the  ultimate  wake  pitch  which 
are  about  89  percent  of  the  measured  values,  corrected  to 
open  water.  The  pitch  predicted  from  equilibrium  of  the 
ultimate  wake  (8)  overestimates  the  measured  value  of  around 
three  percent  for  an  assumed  core  radius  of  1/50  of  the  propeller 
radius. 

The  sensitivity  of  the  calculated  result  to  core  radius 
is  seen  to  be  very  weak.  An  assumed  core  radius  of  1/10 
decreases  the  predicted  pitch  by  three  percent,  thus  coinciding 
almost  exactly  with  the  measured  value.  However,  one  can  find 
no  physical  justification  for  such  a large  core  radius,  except 
for  the  possible  argument  that  this  allows  for  the  lack  of 
complete  roll  up  of  the  initial  vortex  sheet  into  tjhe  core. 

At  the  design  advance  coefficient,  the  actuator  disc 
formula  again  under  predicts  the  pitch,  although  the  results 
are  about  95  percent  of  the  measured  values.  The  pitch  predicted 
from  the  ultimate  wake  model  is  again  slightly  high. 

The  results  for  the  stock  motorboat  propeller  are  given 
in  Table  9.  The  trends  are  remarkably  similar  to  those  of 
Table  8,  with  the  actuator  disc  approximation  slightly  under- 
predicting the  measured  pitch,  and  the  ultimate  wake  model 
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slightly  over-predicting  its  value.  The  fact  that  reasonable 
results  are  obtained  at  high  values  of  the  propeller  thrust 
coefficient  is  encouraging. 


Open  Water  J = 0.6, 
Tunnel  J = 0.648 

MODEL 

NUMRFRC 

7rAi  *=  .960 

4118 

4132 

4133 

4143 

Measured  CT  in  Tunnel 

1.82 

1.60 

2.09 

1.67 

Measured  Pw/D  in  Tunnel 

.850 

.828 

.870 

.791 

Measured  Pw/D  - corrected  to  open 
water 

.918 

.894 

.940 

.854 

Calculated  Pw/D  - Actuator  Disc 

.804 

.784 

.828 

.790 

Calculated  Pw/D  - Ultimate  Wake 
rc  = 1/50 

.949 

Calculated  Pw/D  - r = 1/25 
c 

.933 

Calculated  Pw/D  - r = 1/10 
c 

.912 

Open  Water  J = 0.833 
(DESIGN) 

Tunnel  J = 0.862 
irXi  = 1.041 


Measured  CT  in  Tunnel 

Measured  Pw/D  in  Tunnel 

Measured  Pw/D  - corrected  to  open 
water 

Calculated  Pw/D  - Actuator  Disc 
Calculated  Pw/D  - Ultimate  Wake 


Calculated  Pw/D  - r£  = 1/25 
Calculated  Pw/D  - r£  = 1/10 


.939 

1.029 

1.019 

1.006 


MODEL  NUMBERS 


4132 

4133 

4143 

.561 

.606 

! 

.495 

.963 

.960 

.933 

.996 

.993 

.965 

.937 

.944 

.926 

Table  8.  Summary  of  measured  and  calculated  results  for  DTNSRDC  3-blade  propeller 
series . 


i 


1 

OPEN 

WATER  J - 

1 

i 

.348 

.438 

.544 

.577 

1 i t 

' Measured  C^,  in  Tunnel 

3.650 

1.815 

.824 

! .640 

| 

i 

i Measured  Pw/D  in  Tunnel 

.558 

.584 

! .655 

.651 

' 

1 Measured  Pw/D  - corrected  to  open 

water 

.632 

.633 

.685 

.673 

Calculated  Pw/D  - Actuator  Disc 

.549 

.587 

.639 

.658 

i , Lifting  Line  ttX^ 

/*  i i * j a.  /n  Ultimate  Wake 
Calculated  Pw/D  - _ . 

rc 

.625 

.666 

.651 

.687 

i 

i 

.701 

| 

i 

\ Measured  r /R  | 

1 . « | 

’ 1 

.821 

.803 

.803 

.830 

i Table  9.  Summary  of  measured  and 

calculated  results 

for  12  x 8 

stock  propeller 

VI.  APPLICATION  TO  PROPELLER  LIFTING-SURFACE  DESIGN 

The  propeller  blade  section  design  computer  program 
designated  MIT-PBD-9  [7]  was  written  with  the  provision  for 
modification  of  the  trajectory  of  the  wake  vortices.  The 
wake  model  developed  in  the  present  work  has  been  incorporated 
in  a version  of  this  program  designated  PDB-A. 

The  transition  wake  is  defined  by  a set  of  discrete 
vortices  originating  at  each  grid  point  on  the  trailing  edge 
proceeding  to  specified  hub  and  tip  "roll-up  points".  The 
radius,  as  a function  of  angular  position,  is  simply  assumed  to 
be  two  parabolic  segments  with  zero  slope  both  at  the  trailing 
edge  and  at  the  roll-up  point.  The  pitch  of  each  element  is 
uniquely  determined  from  the  differences  in  angular  and  axial 
coordinates  of  the  trailing  edge  starting  point  and  the 
roll-up  point. 

Downstream  of  the  roll-up  point,  the  system  of  discrete 
vortices  degenerates  to  a single  helical  vortex  from  each  blade 
and  an  axial  hub  vortex.  This,  of  course,  is  simpler  from  a 
computational  point  of  view  than  the  usual  numerical  represent- 
ation of  a complete  helical  vortex  sheet. 

A hypothetical  propeller  with  a blade  outline  identical 
to  4118,  but  with  a design  of  J of  0.6  and  a design  of 

1.816  was  designed  using  both  PBD-9  and  PBD-A,  and  the  results 
are  given  in  Table  10. 


PITCH  RATIO 
P/D 


CAMBER  RATIO 
fo/c 


-I 


r/R 

! H WAKE 
(PBD  - 9) 

DEFORMED 
WAKE  (PBD-A) 

U WAKE 
(PBD  - 9) 

DEFORMED 
WAKE  (PBD-A)  | 

] 

.3 

.912 

.949 

.0240 

.0259 

.5 

.990 

1.008 

.0347 

.0348 

.7 

1.011 

1.009 

.0345 

.0342 

.9 

.996 

.578 

.0354 

.0352 

.95 

.996 

.976 

.0431 

.0428 

The  results  labelled  "gi  WAKE"  correspond  to  a helical 
wake  with  a constant  pitch /diameter  ratio  of  0.960  determined 
by  a Lerbs  induction  factor  method  lifting-line  computer 
program.  The  results  labelled  "DEFORMED  WAKE"  correspond  to 
a rolled-up  wake  with  a pitch/diameter  ratio  of  0.849  and 
a radius  of  83%  of  the  propeller  radius.  The  roll-up  point 
was  specified  to  have  an  axial  coordinate  of  41%  of  the 
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propeller  radius  and  an  angular  coordinate  of  90  from  the 
blade  tip. 

The  tabulated  results  indicate  that  the  influence  of 
wake  deformation  on  blade  shape,  for  a specified  circulation 
distribution,  is  small.  Furthermore,  the  efficiencies  for 
the  two  cases  computed  by  integration  of  forces  on  the  blade 
in  accordance  with  the  procedure  outlined  in  [7]  were 
essentially  identical. 

This  would  indicate  that  there  is  no  particular  reason 
to  depart  from  the  classical  moderately  loaded  wake  model  for 
the  purpose  of  lifting-surface  design  computations. 
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